Do the cosmological observational data prefer phantom dark energy? 
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The dynamics of expansion and large scale structure formation of the Universe are analyzed for 
models with dark energy in the form of a phantom scalar field which initially mimics a A-term and 
evolves slowly to the Big Rip singularity. The discussed model of dark energy has three parameters - 
the density and the equation of state parameter at the current epoch, Q^e and wo, and the asymptotic 
value of the equation of state parameter at a — > 00, c 2 a . Their best-fit values are determined jointly 
with all other cosmological parameters by the MCMC method using observational data on CMB 
anisotropics and polarization, SNe la luminosity distances, BAO measurements and more. Similar 
computations are carried out for ACDM and a quintessence scalar field model of dark energy. It is 
shown that the current data slightly prefer the phantom model, but the differences in the maximum 
likelihoods are not statistically significant. It is also shown that the phantom dark energy with 
monotonically increasing density in future will cause the decay of large scale linear matter density 
perturbations due to the gravitational domination of dark energy perturbations long before the Big 
Rip singularity. 
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I. INTRODUCTION 



Among the large number of dark energy models the 
phantom model has probably the most exotic behavior 
because it violates the null energy condition p + p > 0. 
It was proposed first by R. Caldwell QJ and A. Starobin- 
sky independently in 1999 as a possible explanation 
of the accelerated expansion of the Universe discovered a 
year earlier by two teams 0,0] , measuring the luminosity 
distance-redshift relation from SNe la in distant galaxies. 
It has been shown that this model of dark energy does not 
contradict the cosmological tests based on present data. 
The model was supported later by reconstruction of the 
equation of state parameter using new SN la data § and 
a fully consistent analysis of CMB, large scale structure 
and SN la data [||. The 'doomsday scenario' described 
in Q, caused by the Big Rip singularity which is pre- 
dicted by these models, surpasses the fantasy of science 
fiction writers. Since then, many authors have analyzed 
various aspects of phantom dark energy and confirmed 



its validity as candidate for dark energy [§[. Moreover, 
WMAP [g] data combined with either SNe la or BAO 
(Baryonic Acoustic Oscillations) both prefer the phan- 
tom model of dark energy. This is why such models are 
the subject of active research in recent years: about a 
thousand papers devoted to phantom dark energy can 
be found in the publication databases. 

In Ref. [l| it has been shown that phantom dark energy 
can be modeled by a minimally coupled scalar field with 
classical Lagrangian apart from the kinetic term which 
has the opposite sign. Other scalar field Lagrangians 
with non-canonical kinetic term leading to phantom-like 
properties of dark energy have been considered by differ- 
ent authors later. It turns out that some phantom models 
emerge effectively from the gravity sector of brane- world 
models [Io| . from superstring theory (Tlj . from Brans- 
Dicke scalar-tensor gravity [l2l . [l3| and from quantum 
effects leading to violations of the weak energy condition 
on cosmological scales [3 llEj . Some of these models 
have phantom properties only at the current stage of the 
evolution of the Universe but did not have them at early 
time or they lose this feature in the future 1 . 
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1 We do not consider here models crossing the phantom divide like 
e.g. 'quintom' dark energy, see Ref. [id . 
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It was shown by Carrol et al. (2003) and Cline et 
al. (2004) [8[ that minimally coupled scalar fields with 
a linear negative kinetic term may cause a UV quan- 
tum instability of the vacuum manifesting itself in the 
production of pairs of ghosts, photons or gravitons as 
a consequence of the violation of the null energy condi- 
tion . For late type phantom scalar fields the produced 
ghosts typically carry low energy, so, their decay rates 
are strongly time-dilated. On the other hand, the time 
scale of this instability for phantom dark energy can be 
much larger than the cosmological one, making this effect 
unsuitable for constraining the parameters of the model 
at the present level of observations. This is why in this 
paper we concentrate our attention on the classical prop- 
erties of phantom scalar field models of dark energy and 
on possibilities to determine their parameters by compar- 
ison of predictions with available observational data. 

The main feature of phantom dark energy is its 
strongly negative equation of state parameter, Wde < — lj 
which implies an energy density increasing with time. If 
Wde = const, the energy density of such a field is zero 
at the Big Bang, it starts from "nothing", that is why 
it is dubbed ghost or phantom. In this paper we will 
show that this true phantom can arise as a special case 
of a dynamical scalar field with barotropic equation of 
state which we call PSF (for phantom scalar field with 
barotropic equation of state) in the following. In the 
general case it starts as vacuum energy or a cosmologi- 
cal constant with Wde = — 1 and evolves to a lower value 
Wde < — 1 at the current epoch and to a Big Rip sin- 
gularity in the distant future. The PSF model of dark 
energy has three parameters, we determine their best-fit 
values jointly with other cosmological parameters using 
available observational data on CMB anisotropics and 
polarization, SNe la luminosity distances, BAO measure- 
ments and others. We compare the maximum likelihood 
of the studied PSF model with corresponding values for 
a barotropic scalar field with positive kinetic term, which 
we call QSF (for quintessence scalar field with barotropic 
equation of state), and with ACDM for the same data 
sets. 

The paper is organized as follows. In Section|TT]we dis- 
cuss the cosmological dynamics of phantom scalar fields 
with barotropic EoS. In Section HTll we analyze the grav- 
itational instability of PSF and its effects on structure 
formation. In Section [IVI we present the results of an es- 
timation of PSF parameters and compare the goodness 
of fit of three types of models (PSF , QSF and ACDM) 
for the same data sets. The conclusions are found in 
Section M 



2 This can be prevented by introducing the squared kinetic term 
in the Lagrangian as in the ghost condensate model [ljj or by 
second derivatives of the scalar field as in the kinetic braiding 
scalar-tensor model llql . 



II. EVOLUTION OF PHANTOM SCALAR 
FIELD AND EXPANSION OF THE UNIVERSE 

Let us analyze phantom dark energy from a single min- 
imally coupled scalar field. Real scalar fields with clas- 
sical Lagrangian L — X — U{4>), where X = (p.itfi'' 1 /2 
is the kinetic term, or with a tachyonic Lagrangian, 
L = —U(<f>)y/i — 2X, cannot be models of phantom dark 
energy since the kinetic term is positive. The simplest 
Lagrangian allowing w < — 1 is that of a classical scalar 
field with a kinetic term of opposite sign: 



L„ 



-x-u(<j>). 



(i) 



This has been proposed in [l[ and @ and it is inherent 
from spacelike brane constructions in string theory [L9| . 
In this case the energy density and pressure are the fol- 
lowing linear combinations of X and U: 



Pde = -X + U(<j>), Pde = -X- U(<j>). 



The EoS parameter 



Wde 



-X-U 

-x + u 



(2) 



(3) 



for positive values of X and U is less than — 1. In order 
to explain the accelerated expansion of the Universe at 
the current epoch (q Q < 0) the phantom scalar field must 
satisfy two conditions: 



a) < < U<®, b) U<® + 2X<® > p®/2. 



(4) 



We specify the scalar field model of dark energy by 
the condition Pde/f>de = c\ = const. This is equiva- 
lent of the generalized linear barotropic equation of state 
Pde = c\pde + C, where C is a constant. Represent- 
ing the equation of state as pde — WdePde from the 
energy-momentum conservation law in the Friedmann- 
Robertson- Walker metric with the scale factor a(t) the 
analytic solutions for Wdeifl) and Pdeifl) have been ob- 
tained in 201: 



Wde(a 



Pde 



(l + c 2 a )(l+w ) 



1. 



~ f + w - {wq - c 2 a )a 3( - 1+c ^ 

( o ) (l + w )a-^+^+cl-wo 
Pde l + d 



(5) 
(6) 



where Wq = Wd e (a = 1) is the initial condition for in- 
tegration of the differential equation w' de = 3a _1 (l + 
Wde)(wde — c 2 a ) (a prime denotes the derivative with re- 
spect to the scale factor a). This simplifies significantly 
the analysis of cosmological consequences of this field. 
From ([5]) it follows that at the current epoch (a = 1) 
Wde — wq. For a phantom scalar field (taking into ac- 
count that ap' de = — 3/9^(1 + Wde) > 0) the value of wq 
must be less than -1. The additional condition Wd e < — 1 
for any a gives immediately the constraint c\ < w§ and 
Wde{oo) = c~. If Cq > wq the phantom divide (wde = — 1) 
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FIG. 1: Left column: top panel - the dependence of EoS parameter on the scale factor for barotropic phantom scalar field 
with wo=-l-2 and different c\ (-2.0, -1.3, -1.2, -1.1); bottom panel - the dependence of the dark energy density (in the units of 
critical one at the current epoch) on the scale factor with the same EoS parameters. Right column: top panel - the dynamics 
of the expansion of the Universe with barotropic phantom scalar field: H 2 (a) (top panel) and q(a) (bottom panel) for the same 
models as in the left panels. 



is crossed in the past, the energy density of scalar field 
changes sign and the equation of state parameter diverges 

at a s given by a^ 1+c °' ) = —(1 + ivq)/(c^ — u>o). Here we 
exclude such models from our considerations. In Fig. [T] 
(left column) the evolution of Wd e and pde is shown for 
some values of wq and <r a . One can see that the energy 
density of phantom barotropic scalar field increases for 
all combinations of wq < — 1 and c a < — 1 . 

The quintessence sector of barotropic scalar fields 
(-1 < w Q < -1/3, -1 < c\ < 0), denoted QSF, where 
the density of the scalar field decreases with expansion 
more slowly than the density of matter, has been stud- 
ied in detail in our previous papers |20l - [23 |. This case 
requires c 2 > wq in order for the dark energy to remain 
positive. A comparison of the models is given in Table HI 



TABLE I: The asymptotic behavior of quintessence (QSF) 
and phantom (PSF) barotropic scalar field models. 



QSF 

W d e 
Pde 



PSF 

Wde 
Pde 







CO 



Pd 



-1 

(0) 



a = 1 



Wo 

o {0) 

Pde 



Wo 
Pde 



o (0) - 

Pd, 



1+c 



OC 



Future of the Universe 



-1 < w < -1/3 
w < cl <0 : 

eternal inflation 
-1 < c 2 a < w < -1/3 : 
Big Crunch singularity 



c a < wo < -1 
Big Rip singularity 



We consider a multicomponent model of the Universe 
filled with non-relativistic particles (cold dark matter and 
baryons), relativistic particles (thermal electromagnetic 
radiation and massless neutrino) and a phantom scalar 
field as described above, termed PSF+CDM. The back- 
ground Universe is assumed to be spatially flat, homoge- 
neous and isotropic with Friedmann-Robertson- Walker 
(FRW) metric, 

ds 2 = gijdx l dx-' = a 2 (r])(dr) 2 — 8 a pdx a dx^)^ 

where rj is the conformal time defined by dt — a(ri)dr] and 
01(77) is the scale factor, normalized to 1 at the current 
epoch (here and below we put c = 1). Latin indices i, 
j,... run from to 3 and the Greek ones are used for the 



spatial part of the metric: a, /?, .. = 1, 2, 3. The dynamics 
of the expansion of the Universe can be deduced from the 
Einstein equations 

Rii - \gaR = SttG +T^+ T^) , (7) 

where Ra is the Ricci tensor and t/™^ , T^' , are the 
energy-momentum tensors of non-relativistic matter (m), 
relativistic matter (r), and dark energy (de) respectively. 
Assuming that the interaction between these components 
is only gravitational, each of them satisfies the differential 
energy-momentum conservation law separately: = 
0. Einstein's equations together with energy conservation 
lead to the Friedman equations, which describe the rate 
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FIG. 2: Right panels: The dependences of potentials U (solid line), field variable <f> — (f>i (dotted) and the kinetic term X 
(dashed) as functions of the scale factor a for PSF with Lagrangian JT]) and barotropic EoS (O with different relations between 
wo and c 2 , (c 2 , = —2.0, —1.2, —1.1 from top to bottom, wo = —1.2 for all). Left panel: The dependences of potentials U 
on (<f> — (j>i) for the same PSF models. In all except the bottom right panels the dependences are in log — log scales, in the 
right bottom panel we uselog-lin scales. The potential and the kinetic term are in the units of current critical energy density, 
3c 2 Hq /StvG, the field variable is in units of ^/3c 2 /87rG. In the right panels field evolves from left to right. 



and acceleration of the expansion of the Universe: 



H 

q 



Ho^flr/a* + r> m / a 3 + fl de f(a), (8) 
1 2n r /a 4 + fl m /a 3 + (1 + 3w de )Q de f(a) 



n r /a 4 + Vl m /a 3 + ttdef(a) 



(9) 



where f(a) = pde{a) / Pde(l) ■ Here H = a /a 2 is the Hub- 
ble parameter (expansion rate), q = — (aa/d 2 — l) is the 
deceleration parameter and an overdot denotes deriva- 
tive w.r.t. conformal time rj. Eqs. ((8|)-((9]) completely de- 
scribe the dynamics of expansion of the homogeneous and 
isotropic Universe. In the past it was dominated by ra- 
diation and matter, in the distant future it is dominated 
by the phantom scalar field, as shown in the Fig. [TJ One 
can see that the rate of expansion H (top right panel) 
decreases when the deceleration parameter is q > and 
starts to increase when the phantom component begins 
to dominate. In the case of a quintessence scalar field 
with barotropic EoS H decreases always, asymptotically 
approaching a constant value in the case of c a > wq lead- 
ing to late inflation [2pj . 

Assuming a barotropic equation of state the field vari- 
able, the potential and the kinetic term can be obtained 
in terms of the scale factor and the three model param- 



eters (/0^,uiq,c 2 ) as follows: 



0(a) 



Jo a ( 2 + 2 C "> H (a' ) 
JT( s _ (l~cl)(l + W o)a-^+^+2(cl-w a ) (0 ) 



X 



2(1 + 4) 
(a) = -l±^a-3(^)p(°). 



(10) 



The phantom barotropic scalar field (wq < —1, c 2 , < — 1) 
has real field values if the current dark energy density 
is non-negative. Its kinetic term X(a) is positive for all 
a, the potential U(a) is positive as long as c a < wo- 
lf wo < c a < —1, U(a) starts from the negative value 
(c 2 , — wo)p d ° ! ! /(I + c 2 ) at a = 0, changes the sign at 

a (p=0) = [2{w - c 2 )/(l - c 2 )(l + w )f^fcl), which 
for the phantom case is always < 1. In any case U(a) 
increases with a. This distinguishes the phantom scalar 
field from the quintessence field. The evolution of the 
potential, U(a), field variable 3 , (ft(a), and kinetic term, 
X(a), for models with c 2 < Wq < — 1, c 2 = wq < — 1 and 



3 For definiteness and graphic representation we put the unde- 
tectable initial value <pi equal to zero. 
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the scale factor: 
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FIG. 3: Dependences of the scale factor on time, a(t), for 
cosmological models with quintessence/phantom scalar fields 
with wo = -1 ± 0.2 and e£ = -1 ± 1 (dotted line), -I ± 0.3 
(dashed one), — 1 ± 0.2 (dash-dotted), — 1 ± 0.1 (dash-three- 
dotted). The upper sign is for QSF, the lower one for PSF. For 
the ACDM model a(t) is shown by thick solid line. Note also 
that the limiting behavior for PSF with 1 + <? a /*■ does tend 
to the cosmological constant behavior while 1 + c a /* does 
not. In all models fl m — 0.3, fide = 0.7, Ho = 70 km/s-Mpc. 



wo < (? a < — 1 arc shown in Fig. [5] In the right panels 
the reconstructed potential U(4>— 4>i) is also shown. The 
accelerated expansion of the Universe is caused by the 
rolling of the field up to the maximum of its potential, 
inversely to the case of quintessence scalar field. The 
energy density and pressure are smooth monotonic func- 
tions of a for all relations between c a < — 1 and Wq < — 1 , 
while Wde has a discontinuity of the second kind in the 
case wq < (? a < — 1, when the scalar field energy density 
passes through zero (see Fig. [TJ left bottom panel). 

Another difference between PSF and QSF is in their 
asymptotic behavior: PSF mimics a cosmological con- 
stant at the Big Bang for any < — 1 (wde tends to -1 
when a tends to 0), while QSF does this at a — > oo in 
the case of decreasing EoS parameter, see Tabic fl] PSF 
always starts as cosmological constant with pdeifl = 0) = 

Pd°e( c a ~ w o)/(l + c 1), which is positive for c 2 a < wq and 
negative if Wq < c a < — 1. This property distinguishes 
the barotropic phantom scalar field from the "standard" 
phantom dark energy discussed in 0, 0] , where the den- 
sity starts from zero at a = 0. For PSF only the special 
case with c a = wg has this behavior. 

In the future, when a ^ 1, the energy density of PSF 
increases as pde(a) oc (1 + w ) / '(1 + c^) p d °J a~ 3(1+c ») while 
Wde tends to c a (eqs. (|5|) and ([6])). The repulsion proper- 
ties of PSF increase and in finite time they reach and out- 
match first the gravitational force, then electromagnetic 
forces and finally strong interactions. All bound struc- 
tures in the Universe - galaxies, stars, planets, atoms and 
protons - will be ripped apart in finite time. This singu- 
larity is dubbed the Big Rip Q and the moment when it 
happens can be estimated from the time dependence of 



t{a) 



da 1 



(11) 



which can be computed numerically using ([8]) for any 
cosmological model and parameters of scalar field with 
barotropic EoS. In Fig. [3] we present the time depen- 
dences of scale factors, a(t), for cosmological models with 
PSF with the same parameters as in Fig. [TJ For compar- 
ison we show also a(t) for the ACDM and QSF+CDM 
models with symmetrical values of wq and (? a relative to 
the phantom divide line. Phantom regime for a(t) lies 
above the a(i)-curve for ACDM with the same cosmolog- 
ical parameters, while quintessence range is below. At 
a ^> 1, when radiation and matter terms in (JSJ) can be 
neglected, we obtain the following approximate analytic 
expression for a(t): 



a(t) 



l + ci 



(t - t ) + 1 



(12) 



This shows that a — >• oo is reached within finite time 



tBR 




(1 + W )Q, d e ' 



(13) 



which is noted as the time of Big Rip 0, 0] • This super- 
fast expansion leads also to freezing of the particle hori- 
zon r p at some r™ ax and to a decay of the event horizon 
r e to zero at t tBR- (We denote expansion which is 
faster than exponential by 'super- fast' expansion.) More 
precisely, in comoving coordinates these quantities be- 
have as follows: 



r P (t) = 



da' 



da' 



a' 2 H(a') 



Starting from af 3> 1, when the matter component in ([5]) 
can be neglected, the integral for r p from to a > a/ 
can be computed analytically, 



r„(t) 
I(a f ,a) 



r p(tf) + T (af,a) 

2 r~ 



(1 + 3cl)H Q V (1 + w )Q de V f 



I{af, a) tends to when a/ — > oo. The finite time within 
which the singularity is reached, the freezing of the par- 
ticle horizon and the decay of the event horizon suggest 
that the dynamics of expansion of the Universe domi- 
nated by phantom dark energy is like the free fall into a 
Schwarzschild black hole. 

The positive energy density of PSF becomes infinite 
in finite time f) 13[) . overcoming all other forms of mat- 
ter. The phantom scalar field dark energy rips apart first 
clusters of galaxies, later the Milky Way and other galax- 
ies, then the solar system, a bit later the Sun and stars, 
the Earth and finally "the molecules, atoms, nuclei and 
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nucleons, which we are composed of, before the death of 
the Universe in a Big Rip" (see Table 1 in ]7]). Will this 
be the end of Everything? Maybe this will be the begin- 
ning of new worlds - if PSF reaches the Planck density, 
quantum fluctuations or interaction of the field with par- 
ticles (the phenomenon of confinement) will lead to the 
inflation in some regions at Planck scales. In Ref. 12] it 
has been demonstrated that in a phantom Big Rip quan- 
tum gravity effects might drastically change the future of 
our Universe, removing the singularity in a quite natural 
way. 



III. GRAVITATION INSTABILITY OF THE PSF 
AND LARGE SCALE STRUCTURE FORMATION 

Before determining the PSF parameters let us discuss 
briefly the gravitational instability of such a scalar field 
and its impact on matter clustering. The complete sys- 
tem of evolution equations for cosmological perturbations 
of cold dark matter, baryons, massless and massive neu- 
trinos as well as radiation based on Einstein's equation, 
conservation laws and Boltzmann equations is presented 
in [23|, [24[ • Accurate line of sight integration can be per- 
formed using publicly available codes like CAMB [25| and 
CLASS [11]. Here we use CAMB modified to include the 
expressions ([5]), (O, ([8]) presented above and evolution 
equations for scalar field perturbations [2(|. In the gen- 
eral case the scalar field pressure perturbation includes 
in addition to the adiabatic component 8p^ ad ^ also a non- 
adiabatic mode 5p( nad \ Sp( de ) = c 2 Sp^ d e) + Sp^f\ which 
is interpreted as the manifestation of intrinsic entropy of 
the scalar field (see 20] and references therein). In the 
dark energy rest-frame the total pressure perturbation is 
expressed as 5pd e — c 2 8p de) where the effective sound 
speed 4 for the scalar field with given Lagrangian can 
be calculated as cl = p, x /p,x = L,x I (2XL, X x +L, X ), 
first proposed in [27| • For the Lagrangian (fT|) c 2 = I and 
the effective speed of sound (propagation of the pertur- 
bations) is equal to the speed of light. To understand the 
gravitational instability of PSF and its impact on large 
scale structure formation in the matter and dark energy 
dominated epochs it is sufficient to analyze the following 
set of differential equations: 

Sde + 3(c 2 s - w d e)aH5 de + (1 + w de )^ 



+(l + w de ) 



k + 9a z H 



2 Tj2 



c 2 -c 2 



V de + aH(l - 3c 2 )V de 



k 

c 2 k 



1 + w dt 



Vde = 0, (14) 

■5 de = 0, (15) 



(16) 



h + -h = -8%Ga 2 (p m S m + (1 + 3w de )pdeSde), (17) 

describing the evolution of density perturbations of 
dark energy 5 de = 5p de /p de and matter S m = (5p b + 
Spcdm) I (pb + Pcdm), velocity perturbation V de of dark en- 
ergy and the evolution of metric perturbations h = h\ in 
the synchronous co-moving to dark matter gauge. Here 
as above an overdot denotes the derivative w.r.t. confor- 
mal time 77. Adiabatic initial conditions for matter com- 
ponents can be found in [23| and the initial conditions for 
the early time, when the scalar field is subdominant [20j 
are as follows 



J de 



(4 - 3c 2 s )(l + w de ) 



8 + 6c 2 - Ylwde + 9c 2 (w de - c 2 a ) 



v. 



8 {l 

"ill 



(in) 



£ s kT)v, 



8 + 6c 2 -12w de + 9c 2 {w d e-c 2 a ) 



h (m K 



h {in \ (18) 



For positive matter density perturbation 5 (S m > 0) at 
r]i n , the gravitational potential is negative h < and the 
dark energy density perturbation has the opposite sign 
(Sde < 0) for any Wq, c 2 < —1 and c 2 > 0. The ab- 
solute values of their amplitudes increase oc a on super- 
horizon scales, but the density perturbations of the phan- 
tom scalar field change sign and decay after entering the 
horizon at 77 This is shown in Fig. 21 where the 

evolution of Fourier mode k = 0.05 Mpc -1 of density per- 
turbations for dark matter, baryons and phantom scalar 
field is presented for 2 cases: c 2 < wo and c 2 = w - 

Note also, that PSF perturbations in general do not 
obey a barotropic equation of state even if the back- 
ground does. But as we see here, they never become very 
relevant, we shall therefore not stress this any further. 

In the case of c 2 = wq the absolute value of the initial 
amplitude of Sde is larger than for c 2 < wq, but in both 
cases during the structure formation and at present the 
amplitude of Sde is significantly lower than S m . This im- 
plies that the perturbations of minimally coupled scalar 
fields with initial conditions (fT8)l do not significantly af- 
fect structure formation. 

Nevertheless the parameters of a barotropic scalar field 
can be constrained by large scale structure data, since the 
growth rate of matter density perturbations is sufficiently 
sensitive to them. This is illustrated in Fig. [5l where the 
evolution matter density fluctuations, S m (a), is shown 
for models with PSF dark energy. In order to eliminate 
the fc-dependence caused by the baryonic component at 
small scales and emphasize the influence of dark energy 



4 The terms "pressure", "entropy", "effective sound speed" and 
"adiabatic sound speed" of dark energy are used in the literature 
to denote dark energy properties which correspond only formally 
to the corresponding thermodynamical variables. 



5 In the early Universe for superhorizon scales we include also the 
relativistic components, which dominate the energy density. 
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FIG. 4: The evolution of the Fourier amplitude (k — 0.05 
Mpc -1 ) of density perturbations for cold dark matter (dashed 
line), baryonic matter (dotted) and PSF (solid) for <? a < wo 
(top panel) and for c 2 = wq (bottom panel). 
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FIG. 5: The evolution of matter density perturbations from 
the Dark Ages to the present in sCDM, ACDM, QSF+CDM 
(1: w = -0.8, c 2 a = -0.8; 2: w = -0.8, c 2 a = -0.5) and 
PSF+CDM (1: w = -1.2, c 2 a = -1.2; 2: w = -1.2, c 2 a = 
— 1.5) models. Amplitudes are normalized to 0.1 at z = 10 
(a — 0.1). In models with dark energy Q m = 0.3, Qde ~ 0.7. 



one we have normalized the amplitude of matter den- 
sity perturbations to 0.1 at a = 0.1 (free normalization). 
At this time p m /pde ~ 1000, q « 0.5 and the ampli- 
tudes of all Fourier modes evolve essentially equally. For 
comparison the same variables for QSF, ACDM and the 
standard CDM (sCDM) models are also presented. The 
cosmological model with PSF can be distinguished by the 
amplitude of large scale structure inhomogeneities from 
QSF at the 10% level and from ACDM at the level of a 
few percent for < z < 1. 

Let us analyze the evolution of linear density pertur- 
bations in the future. The first issue for clarification is 
the gravitational instability of PSF in the strongly dark 
energy dominated epoch. We have integrated the system 



of differential equations (|14l) - (fl7|) with initial conditions 
(|T%|) up to a = 200, when expansion is already super-fast 
(see Fig. © and p de /p m ~ 10 8 - 10 10 . The results for the 
PSF density perturbations are shown in Fig. [5] (left col- 
umn) in the log-norm scale for different fc-modes (0.0005, 
0.001, 0.0015, 0.002, 0.0025, 0.005, 0.01, 0.05, 0.1 Mpc" 1 ) 
and two expansion rates, which correspond to the models 
with w = —1-2, c 2 = —1.5 (top panel) and c 2 = —1.2 
(bottom panel). Their amplitudes increase slowly and 
the rate depends on the background expansion rate as 
well as on the wave number. In order to visualize the fc- 
dependence of the growth we remove the dependence of 
the initial conditions fc and renormalize the amplitudes at 
a = 0.1 to Sde{k 7 a = 0.1) = 0.1 = S rn (k 7 a = 0.1), so that 
all fc-modes of the PSF density perturbations in Fig. [5] 
have the same amplitudes at a — 0.1. The growth rate is 
larger for small k in the range a — 0.1 — 10 and it is practi- 
cally the same for all modes at a > 10: 5d e oc a~ 3 ^ 1+Ca ^ 2 . 
For the PSF with wo = —1.2 and c 2 a = —1.5 the ampli- 
tude of the k = 0.1 Mpc -1 mode increases from a = 0.1 
to a = 10 by a factor 693, while the amplitude of the 
k = 0.0005 Mpc- 1 mode increases by 10307. For PSF 
with wo = —1.2 and c 2 = —1.2 these numbers are 9 
and 125 respectively. Since the evolution of the ampli- 
tude of the gravitational potential h is driven by the term 
Pm&m + Pde{^ + 3w de )$de (r.h.s. of eq. [TTJ) , shortly after 
a = 1 the perturbations of PSF become important first on 
the largest scales and later also on smaller scales. They 
affect the evolution of matter density perturbations, as is 
shown in the right hand panels of Fig. [6] At scales with 
k > 0.05 Mpc -1 (lines 8, 9 are superimposed in both 
panels) the amplitudes of matter density perturbations 
in the models with PSF increase from a = 1 to a = 10 
only by a factor ~ 1.3 and freeze at this value. On these 
scales the difference between PSF and the ACDM and 
QSF models for the evolution of matter density pertur- 
bations is inappreciable. In the ACDM and QSF models 
all fc-modes evolve similarly to the line 9. However, on 
scales with fc < 0.05 Mpc" 1 the effect of PSF density 
perturbations on the evolution of matter density pertur- 
bations becomes important: the increase of PSF density 
perturbations causes the decay of matter density pertur- 
bations. The larger the scale of perturbation, the earlier 
its amplitude starts to decay 6 . 

Note that this decay of matter density perturbations 
is caused solely by the influence of phantom scalar field 
perturbations, not by the super-fast expansion of the 
background (at a ~ 2 the rates of expansion in the 
PSF+CDM models are close to those in ACDM and 
QSF+CDM, as it can be seen in Figs. [Hand [3]). Ex- 
cluding the effect of perturbations, the amplitudes of all 
fc-modes freeze as it is shown by line 9. This is not the 
beginning of the Big Rip mentioned above, but its analog 



6 In order to visualize this effect in Fig. [6] we have normalized all 
fc-modes of 5 m to 0.1 at a = 0.1. 
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FIG. 6: The evolution of different Fourier amplitudes of PSF (left column) and matter (right column) density perturbations 
from a = 0.1 to a — 200 for models with wo = —1-2, c a = —1.5 (top panels) and wo = —1.2, c 2 a — —1.2 (bottom panels). The 
rest of parameters are the same as in Figs[4] and [5] The different lines correspond to different wave numbers k (in Mpc -1 ) as 
follows: 1 - 0.0005, 2 - 0.001, 3 - 0.0015, 4 - 0.002, 5 - 0.0025, 6 - 0.005, 7 - 0.01, 8 - 0.05, 9 - 0.1 Mpc" 1 . The amplitudes of all 
fc-modes of Sde are normalized to 5de(k, a = 0.1) = S m (k, a = 0.1) = 0.1 at a = 0.1, for all fc-modes. 



TABLE II: The best-fit values for cosmological parame- 
ters of PSF+CDM model and their la limits from the 
extremal values of the N-dimensional distribution deter- 
mined by the MCMC technique from the combined datasets 
WMAP7+HST+BBN+BAO+SN SDSS SALT2 (pi) and 
WMAP7+HST+BBN+BAO+SN SDSS MLCS2k2 (p 2 ). The 
current Hubble parameter Ho is in units kms -1 Mpc -1 . We 
denote the rescaled energy density of a component X by 
u) X = Slxh 2 . 



Parameters pi p2 



tide 


72+ ' 04 


69+ 05 


Wo 


-1.043±°;°f 


-1.002l°;°° 2 


d 


i 1 9 +0.12 

-i.iz_ 50 


1 iq+o.19 

-1.1»_ 42 


10w 6 


n 9O o+0.016 
u - zio -0 .013 


99 q + 0.014 

u.zzo_ 013 


^cdm 


n l-lc+O.Oll 
U.llO_ 01Q 


0-1191°;°?? 


Ho 


70.4+tg 


67.815:2 


n s 


0.961°:°* 


Ofi+0 03 


log(10 10 A s ) 


3.09lg : gg 


q I I +0.08 

J .- LJ --o.n 


7~re.i 


o.os5±g:8S 


0-0861°:°^ 



- log L 3864.86 3859.30 



for linear perturbations. 



IV. THE BEST-FIT PARAMETERS OF PSF 

Let us estimate the best-fit values of parameters of PSF 
similarly to our previous papers [2CH22| devoted to QSF. 
We use the following datasets: 

1. CMB temperature fluctuations and polarization an- 
gular power spectra from the 7-year WMAP obser- 
vations (hereafter WMAP 7) [Mil]; 

2. Baryon acoustic oscillations in the space distribu- 
tion of galaxies from SDSS DR7 (hereafter BAO) 

0; 

3. Hubble constant measurements from HST (here- 
after HST) [H; 

4. Big Bang Nucleosynthesis prio r on baryon abun- 
dance (hereafter BBN) [32l. |33| : 

5. supernovae la luminosity distances from SDSS 
compilation (hereafter SN SDSS) [38j], determined 
using the SALTpJp (hereafter SN SDSS SALT2) 
and MLCS2k2 35] methods of light curve fitting 
(hereafter SN SDSS MLCS2k2). 

In order to find the best-fit value of parame- 
ters of cosmological model with PSF and their con- 
fidence limits we perform the Markov Chain Monte- 
Carlo (MCMC) analysis for two combined datasets: 
WMAP7+HST+BBN+BAO+SN SDSS SALT2 and 
WMAP7+HST+BBN+BAO+SN SDSS MLCS2k2. We 
use the publicly available package CosmoMC [H, H3] 
including code CAMB 25] for the calculation of the 
model predictions. The results of the estimation of 
the PSF parameters jointly with the minimal set of 
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TABLE III: The best-fit values and la confidential ranges 
from the extremal values of N-dimensional distribution for 
dark energy parameters in ACDM and QSF+CDM deter- 
mined by the Markov Chain Monte Carlo technique using 
2 observational datasets: WMAP7+HST+BBN+BAO+SN 
SDSS SALT2 (h, qi) and WMAP7+HST+BBN+BAO+SN 
SDSS MLCS2k2 (1 2 , q 2 ). The current Hubble parameter H 
is in units kms -1 Mpc -1 . 





ACDM 


ACDM 


QSF+CDM QSF+CDM 


Parameters 


li 


1 2 


qi 




fide 
Wo 

d 

H 


n 7 o+0.03 
u - ' °-0.04 

-1 
-1 

70.4i§;» 


70 +0 ' 04 

u - '"-0.04 
-1 

-1 


n -7Q+0.03 
u - '°-0.05 

n QQfi+° 16 
-0 022+ 022 

u -" zz -0.978 
70.2^4.3 


70+ 04 

"• ' u -0.05 

-0.83±g;?f 
-0.88±g;fl 
66.3± 4 ; 3 7 


-logL 


3864.96 


3859.15 


3865.01 


3857.21 



cosmological parameters for the two sets of obser- 
vational data (WMAP7+HST+BBN+BAO+SN SDSS 
SALT2 and WMAP7+HST+BBN+BAO+SN SDSS 
MLCS2k2) are presented in Table HH We mark the 
sets of best-fit parameters by pi and p 2 . Here 
Pi = (fide, w , c 2 a , fl b , {Idm, H , n s , A s , T rei ). Both 
SN SDSS distance moduli datasets prefer values of wq 
slightly lower than -1. In the past, when a — > 0, 
Wde — > — 1- Hence, the PSFs with parameters pi and 
P2 mimic a A-term from the Big Bang up to the current 
epoch, but, due to the instability of the value Wde = — 1, 
even such a small difference changes drastically the future 
fate of the Universe: in ACDM model the Universe as well 
as existing structures (in principle) are time-unlimited, 
while in the PSF+CDM model it reaches the Big Rip 
singularity in finite time, preceded by the destruction 
of the structure from clusters of galaxies to elementary 
particles. More precisely, in the PSF+CDM with param- 
eters pi this happens in ss 152 Gyrs, with p 2 in m 594 
Gyrs. Long before tsR the particle horizon 7 becomes 
r max _ 18710 Mpc in model pi and _ 192 oo in model 

P2, just wl.3 times larger than the current particle hori- 
zon. 

Let us now compare the best-fit values of cosmo- 
logical parameters pi and p 2 established here with 
similar determinations for ACDM and QSF+CDM 
models. Note firstly that the subset of parameters 
(u>b, w c d m , n s, A sl T rei ) practically does not depend on 
the model of dark energy and SN la fitters, since they 
are determined mainly by WMAP7 data. Therefore we 
will compare the parameters, which depend on them, i. 



e., dark energy ones (fide, wo, c a ) an d Hubble parame- 
ter (Ho), and maximum of likelihoods for the same data 
sets. For the same cosmological model but with A and 
QSF instead PSF the best-fit values of (fide, wq, <? a , Ho) 
are presented in the Table Hill SALT2 fitting of SNe la 
light curves prefers models with lower wo than MLCS2k2 
and the best-fit model with PSF has a slightly lower x 2 
(- log L) than ACDM and QSF+CDM. Although the dif- 
ference is not significant, we can say that the dataset 
WMAP7+HST+BBN+BAO+SN SDSS SALT2 slightly 
prefers phantom models of dark energy. The dataset 
WMAP7+HST+BBN+BAO+SN SDSS MLCS2k2, in 
contrary, slightly prefers the quintessence model. SNe 
la distance moduli determined with use of SALT2 fitter 
predict also higher best-fit value of Hubble parameter Hq 
than the distance moduli for the same SNe determined 
with MLCS2k2 fitter do. In the paper Q the differ- 
ences between 2 methods of light curve fitting, SALT2 
and MLCS2k2, are thoroughly analyzed but convincing 
arguments for one or the other are not given. Using them 
for the same SNe la samples, in [38[ it was found that 
the distance moduli determined with SALT2 fitter prefer 
lower best-fit values of Wde = const than those deter- 
mined with MLCS2k2. Our results support this conclu- 
sion. We have also shown in Ref. [22| that the data on 
SNe la from the SDSS compilation with MLCS2k2 fitter 
allow to constrain c\ in the quintessence range while the 
same data with SALT2 fitting do not. 

Let us finally make use of the newer data on SNe la 
distance moduli from 

• SNLS3 compilation (hereafter SNLS3) Q and 

• Union2.1 compilation (hereafter Union2.1) (ioj 



together with data 
glcZ Dark Energy 
4ll. The 
WMAP7+HST- 



7 At the current epoch r® = 14260 Mpc in the model with pi and 
14170 Mpc in the model with p 2 



on BAO from the Wig- 
Survey (hereafter WiggleZ) 
results for the combined datasets 
-BBN+BAO+WiggleZ+SNLS3 and 
WMAP7+HST+BBN+BAO+WiggleZ+Union2.1 are 
presented in Tables IIVI and [V] correspondingly. 

Both these combined datasets prefer phantom fields, 
with the best-fit values of wo and c 2 a lower than in the case 
of SN SDSS. The best-fit values of H for the phantom 
case are also in these cases higher. Nevertheless for these 
datasets the differences in the — logL between PSF, A 
and QSF are still statistically insignificant. 

Note also that the determination of c a is not suffi- 
ciently reliable for all models and datasets except for QSF 
and WMAP7+HST+BBN+BAO+SN SDSS MLCS2k2 
dataset. The constraints on wo are reliable in all cases. 
For more details see Appendix [S] 

The results of determination of cosmological param- 
eters, especially Ho, fide, Wde and c a , presented in the 
Tables HTllIIIl also indicate certain inconsistency or ten- 
sion between fitters SALT2 and MLCS2k2 applied to the 
same SNe la. It was clearly highlighted and analyzed in 
the papers (38l . |42| , but up to now we have not decisive 
arguments for favor of one from them. 
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TABLE IV: The best-fit values and la confidence ranges of 
the N-dimensional distribution for the dark energy parame- 
ters in QSF+CDM, ACDM and PSF+CDM determined by 
the Markov Chain Monte Carlo technique using the dataset 
WMAP7+HST+BBN+BAO+WiggleZ+SNLS3. The current 
Hubble parameter Ho is in units kms -1 Mpc -1 . 



Parameters QSF+CDM 


ACDM 


PSF+CDM 


fide 
WO 
d 

Ho 


n 79+ 04 

u - ' z -0.04 

n qq4+ 014 

n 72+ - 72 

u - ' z -0.28 

70.1±S;_ 


7T fa04 

u - ' °-0.04 
-1 
-1 

70.3±|;5 


7T f0 - 04 

u - ' °-0.04 

i in +01 ° 

-l.XU_ .27 
r-^»_0.33 

n.5±S:_ 


- log L 


3947.00 


3946.75 


3945.98 



TABLE V: The best-fit values and la confidence ranges of the 
N-dimensional distribution for the dark energy parameters 
in QSF+CDM, ACDM and PSF+CDM determined by the 
Markov Chain Monte Carlo technique using the observational 
dataset WMAP7+HST+BBN+BAO+WiggleZ+Union2.1. 
The current Hubble parameter Ho is in units kms -1 Mpc -1 . 



Parameters QSF+CDM 


ACDM 


PSF+CDM 


fide 
WO 
d 

Ho 


n 7 9+0.03 
u -' z -0.04 

-u.y»o_ 005 

U.OO_ 45 

69.7_f; 5 


72+ ' 04 

u - ' z -0.04 
-1 
-1 

69.8_0 


n 70+0. 03 
u - ' °-0.04 

1 iq+0.13 
-X.XO_ 23 

1 c.4+0-54 
-X.O4_ 09 

71 4+ 4 ' 7 


- log L 


3800.89 


3800.76 


3800.48 



Therefore, at the current level of accuracy of cosmolog- 
ical data we cannot clearly distinguish between the na- 
ture of dark energy studied in this work - quintessence 
scalar field, phantom scalar field or simply A. Hopefully 
the data from observations which are currently under way 
or planned will enable us to establish the type of dark en- 
ergy- 



V. CONCLUSION 

We have analyzed the properties of a phantom scalar 
field with barotropic EoS as possible dark energy candi- 
date. We have studied its effect on the past and future 
dynamics of expansion of the Universe as well as on the 
formation of large scale structure. We show that it mim- 
ics a cosmological constant (or vacuum) model of dark 
energy Wde = — 1 at the Big Bang, with an EoS parameter 
decreasing to wq < — 1 at the current epoch and asymp- 
totically approaching c 2 a < u>o in the future. The param- 
eters Wq and cl together with density one fide are free 
parameters of the phantom barotropic scalar field model 



of dark energy which completely define its physical vari- 
ables and properties. The energy density of such a scalar 
field is always positive and increases with the expansion 
of the Universe. The PSF rolls up the potential with 
expansion, the kinetic term increases monotonically but 
remains always smaller than the potential. Only in the 
special case cl = u>o, Wde = const < — 1 the field starts 
with vanishing energy density. The increasing scalar field 
energy density goes in hand with a super-fast expansion 
of the Universe, which implies a freeze in of the particle 
horizon and the decay of the event horizon. Within finite 
time the Universe reaches a Big Rip singularity. 

The best-fit PSF parameters have been deter- 
mined jointly with all relevant cosmological param- 
eters by the MCMC method using the data sets 
WMAP7+HST+BBN+BAO+SN SDSS SALT2 and 
WMAP7+HST+BBN+BAO+SN SDSS MLCS2k2 (Ta- 
ble [II]). From the early epoch up to now the best-fit PSF 
deviates only slightly from ACDM, but the future of the 
Universe differs significantly. In the PSF model with pi 
parameters the Big Rip singularity occurs 152 Gyrs af- 
ter Big Bang, while in the models with P2 parameters it 
occurs in 594 Gyrs. 

The same computations have been carried out 
for ACDM and a quintessence scalar field model 
of dark energy and it was shown that the dataset 
WMAP7+HST+BBN+BAO+SN SDSS SALT2 slightly 
prefers the phantom model of dark energy, while 
WMAP7+HST+BBN+BAO+SN SDSS MLCS2k2 
slightly prefers the quintessence model. However, the 
differences in the maximum likelihoods are statisti- 
cally insignificant. The same conclusions apply to the 
datasets WMAP7+HST+BBN+BAO+WiggleZ+SNLS3 
and WMAP7+HST+BBN+B AO+ WiggleZ+Union2. 1 . 
We hope that more accurate future observations will 
enable us to distinguish between these models. 
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Appendix A: Estimation of the dark energy 
parameters: details 

We have estimated the dark energy parameters fide, wq 
and c\ jointly with the standard cosmological parameters 
performing a series of CosmoMC runs for different scalar 
fields and datasets. Each run had 8 chains converged to 
R- 1< 0.01. 
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The obtained constraints on the dark energy density 
parameter are reliable. 

For wq and (? a the one-dimensional marginal- 
ized posteriors and mean likelihoods are presented 
in the top and middle panels of Fig. [7] for 
the combined datasets WMAP7+BBN+HST+BAO+SN 
SDSS (MLCS2k2 and SALT2) and Fig. E for 
the datasets WMAP7+BBN+HST+BAO+WiggleZ+SN 
SNLS3 and WMAP7+BBN+HST+BAO+WigglcZ+SN 
Union2.1. We see that the mean likelihoods and pos- 
teriors for u> are close to each other for both phan- 
tom and quintessence fields for all considered datasets. 
The shapes of curves for phantom and quintessence look 
like the parts of a single Gaussian cut by prior at the 
phantom divide wq = — 1 (note that the dependences 
in each panel are normalized to 1 at the maximum, 
so the absolute values of posteriors and likelihoods in 
the neighborhood of -1 can be close for quintessence 
and phantom while the normalized values are suffi- 
ciently different). The peak of Gaussian for the dataset 
WMAP7+BBN+HST+BAO+SN SDSS MLCS2k2 cor- 
responds to the quintessence regime of values while for 



all other datasets to the phantom regime. 

For c 2 a the shapes of marginalized posteriors 
and mean likelihoods are significantly different 
for all cases except for QSF with the dataset 
WMAP7+HST+BBN+BAO+SN SDSS MLCS2k2 
(for this case they have the shape of half-Gaussian with 
the center at (? a — —1). Moreover, in all cases with the 
mentioned exception the shapes of mean likelihoods are 
either asymmetric or almost flat, that is, far from a Gaus- 
sian or half-Gaussian shape. Therefore, the constraints 
on value of (? a are not reliable for all cases except for QSF 
combined with the WMAP7+HST+BBN+BAO+SN 
SDSS MLCS2k2 data. 

In the bottom panels of Figs. [7] and [5] the two- 
dimensional mean likelihood distributions in the plane 
(? a — wq are shown. The solid lines present the la and 
2(7 confidence contours. From these plots it follows that 
the contours are not closed but cut by priors on the 
dark energy type (as well as the mean likelihood distribu- 
tion) . To increase the reliability of the constraints a joint 
consideration of both quintessence and phantom fields is 
needed. This will be the topic of a separate paper. 
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WMAP7+HST+BBN+BA0+SN SDSS SALT2. Bottom: the corresponding two-dimensional mean likelihood distributions 
in the plane <? a — wq. Solid lines show the la and 2a confidence contours. 
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FIG. 8: One-dimensional marginalized posteriors (solid lines) and mean likelihoods (dotted lines) for wo (top pan- 
els) and c 2 a (middle panels). From left to right: PSF and WMAP7+HST+BBN+BAO+WiggleZ+SNLS3, QSF 
and WMAP7+HST+BBN+BAO+WiggleZ+SNLS3, PSF and WMAP7+HST+BBN+BAO+WiggleZ+Union2.1, QSF and 
WMAP7+HST+BBN+BAO+WiggleZ+Union2.1. Bottom: the corresponding two-dimensional mean likelihood distributions 
in the plane c„ — wq. Solid lines show the la and 2a confidence contours. 



